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Abstract 

We introduce the Richelot class of superintegrable systems in TV-dimensions whose n < N 
equations of motion coincide with the Abel equations on n — 1 genus hyperellipic curve. The 
corresponding additional integrals of motion are the second order polynomials of momenta and 
multiseparability of the Richelot superintegrable systems is related with classical theory of covers 
of the hyperelliptic curves. 
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1 Introduction 

In classical mechanics, superintegrable systems are characterized by the fact that they possess more 
than TV integrals of motion functionally independent, globally defined in a 27V-dimensional phase 
space. In particular, when the number of integrals is 27V — 1, the systems are said to be maximally 
superintegrable. The dynamics of these systems is particularly interesting: all bounded orbits are 
closed and periodic [5]. The phase space topology is also very rich: it has the structure of a symplectic 
bifoliation, consisting of the usual Liouvillc-Arnold invariant fibration by Lagrangian tori and of a 
(coisotropic) polar foliation [53]. 

The notion of superintegrability possesses an interesting analog in quantum mechanics. Sommer- 
feld and Bohr were the first to notice that systems allowing separation of variables in more than one 
coordinate system may admit additional integrals of motion. Superintegrable systems show accidental 
degeneracy of the energy levels, which can be removed by taking into account the quantum numbers 
associated to the additional integrals of motion, some of their bound state energy levels may be calcu- 
lated algebraically and the corresponding wave functions are expressed in terms of polynomials. One of 
the best examples of this phenomenon is provided by the harmonic oscillator and the Kepler-Coulomb 
problem. A large number of papers have been published on super-integrability in these last years, most 
of them related with second-order integrals of motion (see [31 [HI [IHl [IH dZl dOl dSl [301 [31] for some 
recent results and an extensive list of references). 

A systematic investigations of superintegrable systems have a very long story, which began in 1761 
when Euler proposed construction of the additional algebraic integral for the differential equation 



where / is an arbitrary quartic [TT]. The corresponding superintegrable Stackel systems have been 
classified in [T7] . 

The Abel theorem may be regarded as a generalization of these Euler results. Remind that the 
Abel equations 

Ui(xi)dxi 

play a pivotal role in classical mechanics and that there are two approaches to investigation of the 
Abel equations associated with the Jacobi and Richelot, respectively (see thirtieth lecture in the 
Jacobi book pl2])- In modern mathematics, the first approach or the Abel- Jacobi map is one of the 
main constructions of algebraic geometry which relates an algebraic curve to its Jacobian variety. The 
second approach yields addition theorems theory, moduli theory (modular equations), cryptography 
and so on. 

The aim of this note is to discuss the Richelot construction of addition integrals for the Abel 
equations and construction of the corresponding A^-dimensional superintegrable systems in classical 
mechanics. We treat only classical superintegrable systems here, though the corresponding results for 
the quantum systems follow easily. 

The paper is organized as follows. In Section II, the main Richelot results are briefly reviewed. 
Then we discuss possible application of these results to classification of the superintegrable Stackel 
systems. In Section III, the classification of superintegrable systems separable in orthogonal coordinate 
systems is treated and solved. Some open problems are discussed in the final Section. 



2 The Richelot superintegrable systems 

In this section we use the original Richelot notations [2S] . 

Let y be the algebraic function of x defined by an equation of the form 

$(x,y) = y™ + /i(x)y™-i + • • • + /™(x) = 0, (2.1) 

where /i(x), . . . , fm{x) are rational polynomials in x. According to the Abel theorem a system of the 
p differential equations 

- — dxi + ■ ■ • + - — dxN = 0, I ~ 1,. . . ,p 

dxi dxN 

have additional algebraic integrals if N > p and if ui, • • • , Up being a set of linearly independent abelian 
integrals of the first kind on algebraic curve ()2.1|) . 

The problem of the determining of these integrals consists only in the expression of the fact that 
xi, . . . , x^r constitute a set belonging to a lot of coresidual sets of places, so we have some determinant 
representations for additional integrals. 

For the particular forms of the curve (|2.ip there are some explicit formulae due by Euler [TT], 
Lagrange [51] , Jacobi [13] , Richelot [5S] , Weierstrass [35] and some other [S] [T] [T3| . 

2.1 The Richelot integrals 

Following to Richelot [2S] we will consider hyperelliptic curve 

y2 = fix) = A2„x2" + ^2n-ix2""i + ■ ■ ■ + A,x + Ao (2.2) 
and the following system of n — 1 differential equations 



dxi 
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Let Qk be the values of x at the branch points of the curve ()2.2p and F{x) = (x — xi)(x — X2) • • • (x — x„), 
then in generic case additional integrals of the Abel equations (12. 3p are equal to 



1 



VfW) 

F'(xi) flfe-xi 



1 



F'{xn) afc-x„ 



F'(xi) F'(x„) 
If ^2n = additional integrals of equations (|2.3p look like 



F{ak) 



(2.4) 



A. 



1 



F'(xi) flfc-xi 



+ 



1 



F'(x„) flfe - x„ 



(2.5) 



There are n—\ functionally independent integrals of motion Ck and, of course, their combinations are 
integrals of motion too. 

Using special combinations of Ck we can avoid calculations of the values ak of x at the branch 
points [131 [551 132] • As an example, in his paper Richelot found the following two algebraic integrals 



1 2 



A2„_l(xi H 1- X„) - A2„(xi H 1- x„) 



(2.6) 



and 



V7C 



xfF'(xi) 



x^F'(x„) 



^1 



Xl 



Xl 



The generating function of additional integrals was proposed by Weierstrass [32], see [2] for detail. 



(2.7) 



2.2 Construction of the Richelot superintegrable systems 

Let as apply the Richelot construction to classification of the superintegrable systems in classical 
mechanics. 

Definition 1 The N -dimensional integrable system is the superintegrable Richelot system if n — 1, 
1 < n < N , equations of motion are the Abel-Richelot equations 12.3]) . 

It's easy to get a lot of such superintegrable Richelot systems in framework of the Jacobi separation 
of variables method, see [T71 [30l [31] . 

Let us start with the maximally superintegrable Richelot systems at = n. In this case con- 
struction consists of the one hyperelliptic curve (|2.2p 

^2 = /(A), where /(A) = As^A^" + A2„_iAf + • • • + A + Aq, (2.8) 

and n arbitrary substitutions 

= ^jilj) Aij = Ujiqj)pj, j = 1, . . . , n, (2.9) 

where p and q are canonical variables {pj ,qi} = Sij . 

The n copies of this hyperelliptic curve and these substitutions give us n separated relations 

p2 u2(g^.) = A2nVj{qj)^" + A2n-lVj{qj)f'-^ + ■■■+ AiVj{qj) +Ao, j = 1, . . . ,n, (2.10) 

where 2rt + 1 coefficients A2n , ■ ■ ■ , Aq are linear functions of n integrals of motion Hi , . . . , Hn and 2?! + 1 
parameters ao, . . . , a2n+i- 
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Solving these separated equations with respect to Hk one gets functionally independent integrals 
of motion 



fc = 1, 



(2.11) 



where Uj{qj) are so-called Stackel potentials and S is the Stackel matrix [77] . 

If Hi is the Hamilton function, then coordinates qj(t, ai, . . . , Q!„) are determined from the Jacobi 
equations 

f Sij{qj)dqj 



and 



E 



Stj{qj)dqj 



= r-t, 



(2.12) 



VELi o^kSkjiqj) - Uj{q.j) 



(2.13) 



where t is the time variable conjugated to the Hamilton function Hi. According to Jacobi |T5] these 
equations are another form of the Abel equations (|l.ip and describe inversion of the corresponding 
Abel map. 

In order to use the Richelot results we have to impose some constraints on the entries of the 
Stackel matrix Skj{qj), which give rise to some restrictions on the coefficients Ak [TTl I30j . 

Namely, if we compare n—1 equations (|2.3p and equations (|2.13p at A = x one gets that the Stakel 
matrix in A variables has to be one of the following matrices 



vn-2 



V 1 



\ n— 1 

A"-' 



1 



\n — l 
\n-2 

1 ) 



n, n + 1, 



,2n. 



(2.14) 



(2.15) 



so that 

2n 

= /(A) = A^ffi + A"-iif„_i • • • + Hn-i\ + H,,+Y, A^ . 

j=o 

Such as k is arbitrary number from n to 2n we have a family of the dual Stackel systems associated with 
one hyperelliptic curve l|2.8p and different blocks of the corresponding Brill-Noether matrix [28l [29] 



Remark 1 For any two dual systems with Hamiltonians Hi and Hi the corresponding Stackel matrices 
S''^' and S"'*^) are distinguished on the first row only. These Stackel systems are related by canonical 
transformation of the time t ^ t: 



Hi=y{q)Hi, dt — Y{q) dt, where v{q) 



det 5'('=) 



(2.16) 



det S-e^) 

Such dual systems have common trajectories with different parametrization by the time |29[|20| . Exis- 
tence of the such dual systems is related with the fact that the Abel map is surjective and generically 
injective. 



Remark 2 For the dual systems the corresponding hyperelliptic curves (j2.15p are related by per- 
mutation of one of the a's and Hamiltonian Hi and, therefore, such transformations are called the 
coupling constant metamorphoses |S1 [TSl [55] . Such transformations are related with the reciprocal 
transformations as well [T]. 

Remark 3 There exist the Richelot supcrintegrable systems that can be solved via separation of 
variables in more than one coordinate system. These systems are associated with non-isomorphic 
curves whose Jacobians are isomorphic to one another (either Jacobian of (|2.2p could be isomorphic 
to a strata of another Jacobian or Jacobian of (|2.2p could be isogenous to a product of some different 
curves etc). Such curves have already occurred in the work of Hermite, Goursat, Burkhardt, Brioschi, 
and Bolza, see Krazer \TP and a lot of modern works on the Frey-Kani covers. 
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Now let as briefly consider construction of the superintegrable Richelot systems for which n — 1 
equations of motion among the N equations of motion are the Abel-Richelot equations only. In this 
case to n separated relations (|2.10p have to be complimented by — n separated relations 

^m{Pm,qm,Hi,...,HN) ^0, n< m < N . 

Solving this complete set of the separated equations with respect to integrals of motion Hk we have 
to get TV functionally independent integrals of motion p. lip . As above the Abel equations have to 
coincide with the Richelot equations (|2.3p and, therefore, the nxn block of the N x N Stackel matrix 
has to be matrix as (|2.14[) . If we take into account all these restrictions one gets complete classifications 
of the superintegrable Stackel-Richelot systems. 

The main problem is that we want to get Hamiltonians Hj in some physical variables x instead of 
Hamiltonians (|2.1ip in terms of the abstract separated variables q. According to JJ_, ,30i ^T] it leads 
to some additional restrictions on the coefficients Aj in (|2.8p and substitutions (|2.9p . 

It easy to see that the Stackel integrals of motion Hk (|2.1ip and the Richelot additional integrals 
of motion are the second order polynomials in momenta 

2 

- ^2n-l(vi + • • • + Vn) - A2n{vi + ■■■+ V^f (2.17) 

and 

v\vl ■ ■ - vl - A^i— + ■ ■ ■ ^ ^ - A^i— ^ ■ ■ ■ + ^ . 

(2.18) 

Here Uj and Vj are functions on coordinates only. 

So, in the Stackel-Richelot case all the integrals of motion are the second order polynomials 
in momenta and it allows us to find natural Hamiltonian superintegrable systems on the Riemannian 
manifolds using well-studied theory of the orthogonal coordinate systems and the corresponding Killing 
tensors [1 H [H [H] . 



uipi 
F'ivi) 



^nPn 



F'iVn) 



v^F'iv,) 



^nPn 



3 The Richelot systems separable in orthogonal coordinate 
systems 

All the orthogonal separable coordinate systems can be viewed as an orthogonal sum of certain basic 
coordinate systems [H [9j [HJ [24] . Below we consider some of these basic coordinate systems in the 
n-dimcnsional Euclidean space only. 



3.1 The basic orthogonal coordinate systems 

Definition 2 The elliptic coordinate system {qi} in the N -dimensional Euclidean space Ejv with pa- 
rameters ei < 62 < ■ ■ ■ < cn is defined through the equation 

t^i^-^^ nti(A-e,) 

The defining equation (j3.19p should be interpreted as an identity with respect to A. 

It is possible to degenerate the elliptic coordinate systems in a proper way by letting two or 
more of the parameters coincide. Then the ellipsoid will become a spheroid, or even a sphere if 
all parameters coincide. Rotational symmetry of dimension m is thus introduced if m + 1 parameters 
coincide. 



Example 1 As an example when ei =62, we have 



2 ^2 T[^^^l — 

e(A) = l + -^+V-^= iJ-\ r'=xi+xi. (3.20) 



^-^1 U^-^^ nf=/(A 
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It defines elliptic coordinate system in Eat^i = {r,X3,--- ,xn}- In order to get an orthogonal co- 
ordinate system {qi, • • • , qj^} in E^v, we could complement r with an angular coordinate gjv in the 
{xi, a;2}-plane, for instance through 



Xi=rcosqN, a;2 = r sin g^v , where r = ^resj^j^^^^ e(A) . (3-21) 
At = 3 these equations define the prolate spherical coordinate system. 

When ei = 62 = • ■ ■ = e„ the only remaining coordinate is r = -^/^ xf and — 1 angular coordinates 
have to be introduced on the unit sphere According to [19] these angular coordinates are 

so-called ignorable coordinates. 

Definition 3 The parabolic coordinate system {qi} in E^r with parameters ei < 62 < • • • < ejv-i is 
defined through the equation 

«(A).A-2.„-"f^.5fe|^. (3.22) 

This orthogonal coordinate system can, in fact, be derived from the elliptic coordinate system as well. 
Namely, substitute 



into the ()3.19p and let cn tend to infinity, then drop the primes one gets the parabolic coordinate 
system. 

The parabolic coordinate system can be degenerated in the same way as the elliptic coordinate 
system. 

Example 2 If ei = 62, we have 

As above in order to get an orthogonal coordinate system {gi, • • • , gn} in E^v, we could complement 
r with an angular or ignorable coordinate q^ in the {xi, a;2}-plane defined by p.2ip . At A^ = 3 it is 
so-called rotational parabolic coordinates. 

Definition 4 The elliptic coordinate system {qi\ on the sphere S^r with parameters ei < 62 < • • • < 
cat+i is defined through the equation 



Notice that (|3.24p implies X^S^ — 1- the similar manner we can define elliptic coordinate 
system {qi\ on the hyperboloid Hat with xf, — X^iLi ^^f = 1 [S]- As above these coordinates can be 
degenerated by letting some, but not all, parameters coincide. 

Remark 4 There are some algorithms [H [IJ and software [TB] that for a given natural Hamilton 
function H — T + V determine if separation coordinates exist, and in that case, show how to construct 
them, i.e. how to get determining function e(A). 
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3.2 The maximally superintegrable Richelot systems 

The basic orthogonal coordinate systems is defined by the function 



e(A) = n£iiA^ = ^ Af = 7V,iV±l, 



which is the ratio of the foUowing polynomials 

AT 



U{X) 



(3.25) 



M 



0(A) = n(^-'?j), and «(A)==[](A 

i=i j=i 



(3.26) 



We can describe the maximally superintegrable Richelot systems separable in these coordinate 
systems using the following Proposition. 



Proposition 1 If n ~ N separated relations have the following form 

/ N \ 



a(A) 



4=2 



2N 
J=0 



(3.27) 



where a{\) is arbitrary polynomial, then equations of motion h2.13\) are the Abel-Richelot equations 
If k — n the corresponding maximally superintegrable Hamiltonian Hi 



N 



Hi=T + V =J2 



res 



res 



a(A) 



u2(A)e(A) 



has a natural form in Cartesian coordinates in E„ 

N M 

2 



N M 



i=0 



a(A) 



u2(A)e(A) 



(3.28) 



Here we introduce additional parameter cq — oo. 

If k > n then i^j*^^"-* = v{x) Hi, where function v(x) is defined by h2.16\) . 

It is easy to prove, that these maximally superintegrable Richelot systems coincide with the well-known 
superintegrable systems (U [U [TOl [Ml [20l [26] . 1 For elliptic coordinate system in E^r equation (|3.28p 
yields the following potential 



N 



V = a2N{xl H xl) + ^^, 



x'^ 

For parabolic coordinate system in Ejv one gets 

N-l 



OtiCi) 



i=l 



V = Oi2N{x\ ^ 4x^) + ^NXN + X! IN = ^OL2N ^ 6,, + 2oi2N -\ 



For elliptic coordinate system on the sphere Sat or on the hyperboloid Hat we obtain 

7V+1 



7» 



a(e») 
□^^^(ei-ej)^ 
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Example 3 Let us consider parabolic coordinates (51,52,93) defined by 



e(A) ^ X-2x3 



A — ei A — 62 
whereas the corresponding momenta are equal to 



(A - gi)(A - q2){\ - 53) 
(A-ei)(A-e2) 



Pi 



XlPxi 



X2PX2 



+ 



PX3 



2{q,-ei) 2{q,-e2) ' 2 
In this case the separated relations p.27ll3.3^ look like 

1 



^= 1, 



,3. 



PMi-eiY{qi-e2Y 



[HiX' + H2X + i/3)(A - ei)(A - 62) - a(A) 



\=qi 



,3. (3.29) 



Solving these equations with respect to Hk one gets integral of motion and the following Hamilton 
function 

Tj Pxi Px2 ~^ Px3 / 2 , 2 , A 2\ , I 1 , /„ „„N 

Hi = h ae[Xi + ^2 + 4x3) + 73X3 + — + — + const . (3.30) 

It is maximally superintegrable Hamiltonian with the Stackel integrals of motion H2 , -^3 and two 
additional Richelot integral of motion Ki,2 



(gi - ei)(gi - £2)^1 ^ (g2 - ei)(g2 - 62)^2 ^ (g3 - ei)(g3 - 62)^3 

(gl - g2)(gl - g3) (g2 - gl)(g2 - gs) (g3 - gl)(g3 - g2) 

Yigl +g2 +g3) + Y^^^l +^2 + g3) , 



(3.31) 



K2 



222 

glg2g3 



(gl - ei)(gi - 62)^1 ^ (g2 - ei)(g2 " 62)^2 ^ (g3 - ei)(g3 - 62)^3 

(gl - g2)(gl - g3)g? (g2 - gl)(g2 - g3)gi (g3 - gl)(g3 - g2)gi 

^ g3ei + (g3-g2ei)e2 / 1 , ^ , ^^ 6162^3 /! ^ 1 ^ 1 
2 Vgi g2 g3/ 2 Vgi 92 g3 

In physical variables {x,px) these integrals have more complicated structure. 

It is easy to prove that integrals Hi, H2, H3 and Ki, K2 are functionally independent, of course, 
all these integrals of motion may be obtained in framework of the Weierstrass approach [32] as well. 

Example 4 Now let us consider dual Stackel system and put fc = n + 1 in the Stackel matrix (|2.14p 
from the previous Example. It means that we change one of the coefficients in the separated relations 
(|3.29|) and consider the following separated relations 



Piili - ei)^(gi - 62)^ 



1 r 



(i/iA^ + H2\ + i/3)(A - ei)(A - 62) - a(A) 



Solving these equations one gets superintegrable system with the Hamiltonian 

1 



Hi = y{q) Hi 



2X3 + 61+62 



Hi, 



where Hi is given by (|3.30p . Of course, this canonical transformation of time changes additional 
integrals of motion i^i 2 (|3.30p . 



3.3 The superintegrable Richelot systems 

Now let us consider degenerate coordinate systems for which two or more of the parameters Cj coincide. 

In terms of the separated coordinates defining function e(A) remains meromorphic function with 
n simple roots and m = n, n ± 1 simple poles. For the construction of the Richelot systems we need 
degenerations such that 1 < n < N. 
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In this case in order to get superintegrable Richelot systems with n — 1 additional integrals of 
motion we have to take n separated relations (|3.27p 



p1 u{qif 



i=2 



1 ■r-^ u(X) 
j=n+l °^ ^ ' 



and N — n separated relations for ignorable variables 



j = n + 1, 



,iV. 



(3.32) 



(3.33) 



Here polynomials gj(A) depend on degree of degeneracy and definition of the ignorable variables [HHS], 
whereas Uj{qj) are arbitrary functions on these ignorable (angular) variables qj. 

Solving these equations with respect to integrals of motion Hj one gets the Hamilton function in 
the same form as (j3.28p in which, roughly speaking, trailing coefficient of the polynomial a{X) depends 
on ignorable variables. 

Proposition 2 For degenerate elliptic or parabolic coordinates superintegrable potentials have the fol- 
lowing form \3.28\) 

a{X) - U, 



res 



i=0 



\=ei 



M2(A)e(A)^ 



Co = 00, 



(3.34) 



where Ui — for single roots ei of initial function (A — ei) • • • (A — eu) US. 26\) after degeneration 
Cfc — Cj . For degenerate roots — Cj potential Ui are arbitrary functions on the corresponding 
ignorable variables. 

It allows us classify all the superintegrable Richelot systems using known classification of the orthogonal 
coordinate systems [1 H [TOl H HZl HOI US] ■ 

Example 5 Let us consider prolate spherical coordinate system (91,(72,93) defined by 



e(A) = l + ^±^+ 



^ (A - gi)(A - 92) 
A - ei ' A - 63 (A - ei)(A - 63) ' 

The corresponding momenta are 



qs = arctan I — 



Pi 



XlPxi + X2Px2 



X3PX3 



P2 



XlPxi + X2PX2 



XsPx 



2(91-61) ' 2(91-63)' 2(92-61) 2(92-63)' 

In this case g(A) = (63 — ei)~^(A — ei) and the separated relations (|3. 32113. 33| look like 



P3 = X2PX1 - XlPx2 



-ei)2(9z -63)^ 



P3 = 2(C/(93)-i?3) , 



(HiX + H2){X - ei)(A - 62) - a(A) + 



(A - 63) (63 - ei)iJ3 



where a(A) = a^X'^ + ce^X^ + a2X^ + aiX + uq. 

Solving these equations with respect to Hk one gets integrals of motion and the following Hamilton 
function 



Hi = P^^+P^^^+P^^ + a,{xl +xl+ 4) + 



2x , \X2 



^ - 2a4(e3 + ei) - 03 , 
X3 



where 



71,3 



a(ei^3) 
(ei - 63)^ 
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It is superintegrable Hamiltonian with the Stackel integrals of motion H2 , ffs and additional Richelot 
integral of motion Ki (I2.17p . which is equal to 

^ _ / (gi - ei)(gi - 63)^1 ^ (g2 - ei)(g2 - e3)p2 y (^i - a3)(gi + g2) ^ a4(gi+g2)^ 
V gi - g2 g2 - gi / 2 2 

In physical variables (a;,pa;) one gets the following expression for this integral of motion 

Ki = \ \ai{ei + 63 - - X2 - X3) + as - . 

The second Richelit integral K2 (I2.18P looks like 
K2^( ^ (g2 -ei)(g2 -63)^2 \ ^ ^2^2 _ n ^ _ [^1 _^ 

^ V (gl-g2)g? (g2-gl)gi / Vgl g2/ ° Vgl g2/ 

where 

^1 = i (eie3iJi - (ei + 63)^2 + (ei - 63)^3 - ai) , ^^^ \ (^^^3H2 - 63(61 - 63)^3 - ao) . 

Of course, substituting Hi, . . . , into K2 one gets that Ki = K2, because in this case we have only 
one Abel-Richelot equation, i.e. n — 1 = 1. It means that Hamiltonian Hi in E3 does not maximally 
superintegrable. 

Example 6 Let us consider rotational parabolic coordinates (5i,g2,g3) defined by 

(\\ \ o xj+xj (A-gi)(A-g2) , f xi\ 

6(A) = A - 2a;3 r = r , qs = arctan — , 

A - 61 A -61 \^2/ 

whereas the corresponding momenta look like 

XlPxi + X2PX2 , PX3 XlPxi + X2PX2 , PX3 

Pi = ^TT ^ + ^' = ^ + ^' P3 ^ X2Px, ~ XiPx2 ■ 

2(gi-6i) 2 2(^2 -61) 2 

In this case g(A) = (A — ei) and the separated relations (|3. 32113. 3^ are equal to 



P?.2(gi,2-6i)2 = i 



pI = 2{U{q3) - ^3) 



(i/iA + i72)(A-6i)-a(A) + ^ 



where a{X) = aiX"^ + a3A"^ + Q!2A2 + aiA + ao- 

Solving these equations with respect to Hk one gets integrals of motion and the following Hamiltc 
function 

a{ei)-u(^ 

Hi = + a^ixj +xl + ixj) + 2(2a46i + a3)a;3 + 2 7^ " 3^46? - 2a36i - a2 • 

It is superintegrable Hamiltonian with the Stackel integrals of motion H2 , H^ and additional Richelot 
integral of motion Ki (I2.17p . which is equal to 



(gl - ei)pi ^ (g2 - 6i)p2 
gl - g2 g2 - gl 



j + Y(gl+g2) + Y(gl+g2)' 



Px3 , o 2 , I ^ , 61(0461 +a3) 

+ 204X3 + (2a46i + Q!3)x3 H . (3.35) 



As above Ki = K2 (|2.17ll2.18p in this case. 
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Example 7 Let us consider degenerate elliptic coordinate system on the sphere §3 in E4, so that 
coordinates (91,92, 93) are defined by 



{X~ qi){X~ q2) 



(73 = arctan I — 

X2 



X — ei A — 63 A — 64 (A — ei)(A — e3)(A — 64) ' 

It means that radius of the sphere is equal to i? = J2i=i = 1- 

In this case g(A) = (63 — ei)~^{ei — e4)~^(A — ei) and pair of the separated relations have the 
common form 



(i/iA + H2KX - ei)(A - e3)(A - 64) - a(A) 



+ (63 - ei)(ei - e4)(A - e3)(A - e4)i?3 



(3.36) 



where a (A) is fourth order polynomial with arbitrary coefficients and third separated relation is equal 
to 

pI = 2([/(g3) - H3) . 

Solving separated equations with respect to i/j, one gets integrals of motion and the following Hamilton 
function 



7i + C^(t 



3j ~\~ 



73 , 74 «4 

2;2 ~ i? ' 



7i 



a{ei) 



1-3 ^4 ^1- rij^itilsj ^i)^ 

It is superintegrable Hamiltonian and additional Richelot integrals of motion looks like 



Ki = 



(gi - ei)(9i - e3)(9i - 64)^1 _^ (92 - ei)(92 - e3)(92 - £4)^2 
91 - 92 92 - 91 

(ei + 63 + e4)Hi + - H2 f , \ a^- Hi 



91 + 92 



91 + 92 



(3.37) 



In this case n = 2 and, therefore, Ki = K2 

In this case change of the time (|2.16l) at = ?i + 1 yields the following transformation of pair of 
the separated relation p.32p - (|3.36p 



1 



.(A) • (i/i A^ + H2) - «(A) + i4lTi?3 



In the right hand side of this equations we obtain 2n + 1-order polynomial in A and, therefore, the 
corresponding pair of the Abel equations are no longer the Richelot equations (|2.3p . This change of 
the time preserves integrability, but destroys superintegrability. 



4 Conclusion 

According to [T71 [3D1 [31] there are two classes of superintegrable systems for which the angle variables 
are either logarithmic or elliptic functions. In the both cases one gets additional single- valued integrals 
of motion using addition theorems, which are particular cases of the Abel theorem. 

The main aim of this note is to discuss one of the oldest but almost completely forgotten in mod- 
ern literature Richelot 's approach to construction and to investigation of the superintegrable systems 
separable in orthogonal coordinate systems. Of course, these n-dimensional superintegrable systems 
may be obtained using another known methods (see [3l [H [TOl [Mj [20l [26] and references within). 
Nevertheless we think that new definition (|3.28p . (|3.34p 

(\\ ^'^ 
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of the superintegrable potentials through defining function e(A) of coordinate system and arbitrary 
polynomial a{X) may be useful in applications. 

It will be interesting to get quantum counterparts of the Richelot integrals of motion and to study 
the algebra of integrals of motion in the algebro-geometric terms. Another perspective consists in the 
classification of the Richelot superintegrable systems on the Darboux spaces. 

One more important issue concerns relation of multiseparability of the Richelot superintegrable 
systems with classical theory of covers of the hyperelliptic curves. 
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